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Abstract 

The absence of gravitational radiation in Kinnersley's "photon 
rocket" solution of Einstein's equations is clarified by studying the 
mathematically well-defined problem of point-like photon rockets in 
Minkowski space (i.e. massive particles emitting null fluid anisotro- 
pically and accelerating because of the recoil). We explicitly compute 
the (uniquely defined) linearized retarded gravitational waves emitted 
by such objects, which are the coherent superposition of the gravita- 
tional waves generated by the motion of the massive point-like rocket 
and of those generated by the energy-momentum distribution of the 
photon fluid. In the special case (corresponding to Kinnersley's solu- 
tion) where the anisotropy of the photon emission is purely dipolar we 
find that the gravitational wave amplitude generated by the energy- 
momentum of the photons exactly cancels the usual 1/r gravitational 
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wave amplitude generated by the accelerated motion of the rocket. 
More general photon anisotropies would, however, generate genuine 
gravitational radiation at infinity. Our explicit calculations show the 
compatibility between the non-radiative character of Kinnersley's so- 
lution and the currently used gravitational wave generation formalisms 
based on post-Minkowskian perturbation theory. 
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1 Introduction 



In 1969 Kinnersley constructed an exact solution of Einstein's equations 
containing four arbitrary functions of time. In the literature (see e.g. f^) this 
solution is referred to as a "photon rocket" and is interpreted as a "particle 
emitting null fluid anisotropically, and accelerating because of the recoil" 
Recently, Bonnor pointed out a paradox of this interpretation: one would 
expect the presence of gravitational radiation emitted by the "accelerated 
particle", while, as first shown by Kinnersley [Q], and recently confirmed by 
Bonnor 0, the spacetime described by Kinnersley's solution is asymptoti- 
cally fiat and (though non-stationary) contains no gravitational radiation at 
infinity. Technically, if, say, the "mass function" m(s) is constant for s < Sq 
and becomes smoothly variable for s > sq, Kinnersley's spacetime is identical 
to Schwarzschild's spacetime near past null infinity, and has a (everywhere 
Petrov type D) Weyl tensor decreasing like (in suitable, radiative, co- 
ordinates) near future null infinity, with a Ricci tensor decreasing like 
corresponding to the energy-momentum of the outgoing photon fluid (see 
equations (11) and (12) of 0). 

One possible reaction to this paradox is to remark that the "interpreta- 
tion" on which it is based cannot be given any rigourous meaning within the 
exact, curved spacetime framework of Einstein's theory. Indeed, the notion 
of "point particle" has no mathematical meaning in general relativity, and 
even if we somehow "fill in" Kinnersley's solution to deal with a regular, 
extended massive source, the notion of the "position" of the source (and of 
its "quadrupole moment") will have no unique definition, so that there will 
be no unambiguous way of saying that some massive object is accelerated 
and "must" therefore radiate at future null infinity^. We believe that these 
ambiguities in the interpretation of Kinnersley's solution are indeed very se- 
rious. However, we shall here tackle directly the physics behind the paradox 
raised by Bonnor by noticing that there is a mathematically well-defined si- 
tuation in which we can self-consistently describe a recoiling photon rocket 
and compute the retarded gravitational waves it generates. More precisely, 
we consider, within the framework of special relativity, a point-like massive 
particle (described by a 5 function on a time-like worldline) emitting null 
fluid anisotropically. Working rigourously within the framework of distribu- 
tion theory we prove in section 3 (after having introduced our notation 

^Remembering that Kinnersley's solution can be taken to coincide with a static 
Schwarzschild spacetime before some retarded time, we do not have to worry about the 
possible presence of incoming gravitational radiation. 
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for various quantities associated with a generic worldline in section 2) that 
the total energy-momentum tensor T^,^ = T^^'^'^^^ + T^photons) -g gyery- 
where conserved when the particle loses mass and recoils as expected from 
the fluxes of photon energy and momentum at infinity. In section 4, we then 
compute exphcitly the (retarded) linearized gravitational field generated by 
the total T^jy. We find that the amount of gravitational radiation at infinity 
depends on the anisotropics of multipole order £ > 2 in the photon flux, 
there being a cancellation between the gravitational wave amplitude emitted 
by the energy-momentum distribution of the monopolar and dipolar photon 
flux, and the gravitational wave amplitude emitted by the accelerated mas- 
sive particle. We conclude in section 5 that Bonnor's paradox^ contains some 
truth, but that there is no incompatibility between the non-radiative charac- 
ter of Kinnersley's solution (with its purely dipolar photon flux anisotropy) 
and the consistency of the standard gravitational wave generation formalisms 
in general relativity, based on post-Minkowskian perturbation theory (such 
as a recently developed multipolar-post-Minkowskian formalism [§], 
0' 1131) • Concerns recently raised in the hterature [|1T| about the reliabi- 



lity of the current predictions for the gravitational waves generated by, e.g. 
binary systems are therefore totally unjustified. 



2 Geometry of Minkowski worldlines 

Let C denote a timelike worldline z'^ls), parametrized by the proper time 
s, in Minkowski space (R*^,?]^,^) with 77^,^ = diag(— 1, The four- 
velocity along C is denoted u^{s) = dz^{s)/ds and satisfies rj^^, = — 1. 
To each point in Minkowski space one associates a unique "retarded" 
point Zj^i^x'^) on C, with parameter sr{x^), defined as the intersection of C 
with the past light cone with vertex at x^. In equations, z'^lx) = z^{sji{x)) 
with the function sr(x) defined as the only solution of 

Vf.. (x^ - [sB.ix)]) ix-" - z- \sr{x)\) = (2.1) 

such that x^ — z'ji > 0. One also introduces the null vectors 

y^'ix) = - 4(x) , (2.2a) 

k'^ix) ^ , (2.2b) 

Pr[x) 



^Remember that a paradox is a "statement that seems to say something opposite to 
common sense or the truth; but which may contain a truth" 
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where 

Pr{x) ^ -{x^ - 4) < (2.3) 

denotes the retarded distance between x^^ and the worldhne C Here and in 
the following the Minkowski metric is used to raise or lower indices. Diffe- 
rentiating the defining equation (2.1) yields the spacetime derivative of sr{x), 

d, sr(x) = = -k,{x) , (2.4) 

Pr{x) 

from which follows 

d, y''^S; + y, u'' = S; + k, , (2.5) 

du, PR = p'^ yn-u^ + p^ (yii) Vi^^ k^-Ui^ + pniku) k^ . (2.6) 

Here and in the following ii^ = du^/ ds denotes the proper acceleration of the 
worldline C and we abbreviate Minkowski scalar products as (ah) = a^h^ = 
r]^j^y h" . It is also convenient to introduce notations for the unit spatial 
vector representing the direction of fc'* in the 3- plane orthogonal to m'^, 

k" ^u^ + n" , with (un) = , (nn) = +1 , (2.7) 

and for the projection operator orthogonal to u^^: 

^ixu = Vnu + u^u^ . (2.8) 

With this notation, and the results (2.4)-(2.6), we get also (suppressing the 
"retarded" label R when there is no ambiguity) 

d^i P = rif^ + p{nu) k^ , (2.9a) 

K = p'^i^ixv - n^) - {nil) k^ K , (2.9b) 
n„ = p~^(A^^ - n^) + k^ iii, - {nil) k^ k„ . (2.9c) 
As is well-known the null congruence k^^ is geodesic, shear-free and expanding: 

k" d^k^^Q , (2.10a) 

k^" = -□ sr{x) = - , (2.10b) 
P 

9" k^ k, = p-^[A^u - n,] , (2.10c) 
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k" d, K = 2p-' = ^ {d, . (2.10d) 

Finally, let us note that we can introduce "retarded" polar coordinates cen- 
tered on the worldline C, say {s,p,9,{p), such that the polar coordinates of 
are sr{x), Pr{x) and two angles parametrizing the unit spatial vector 
n^{x), chosen, say, so as to induce the standard metric p'^{d6'^ + sin^6' dip^) 
on the 2-spheres s = const, p = const. This still leaves a large freedom in 
the rotation rate of the basis vectors dn^/d6, (sin^^)^^ dn^^ /dip. We shall not 
need to restrict this freedom. One finds that the Minkowski volume element 
takes (independently of the acceleration of the worldline) the familiar form 

d'^x = p^ ds dp dfl^ , (2-11) 

with dQn = sin 6 dO d(p, while the surface element of the retarded tubes p = 
const, reads 

dSf, = PR p% ds dQn ■ (2-12) 

3 Photon rockets in Minkowski space 

A physical photon rocket would be an extended massive object emitting (in 
the WKB approximation) null fluid anisotropically. Let us prove that one can 
consistently define, in Minkowski space, the point-like limit of such an object. 
We consider an energy-momentum distribution which is partly concentrated 
as a (5-distribution on a time-like worldline C, and partly distributed as a 
null fluid spurting out of C Say 

T''ix)=TlZ^{x)+Ti;:^{x) , (3.1) 

with a "matter" part 

TlZ){x) = Jds m{s) u^{s) u^is) 6^{x^ - z^{s)) , (3.2) 
and a "photon" part 

It has been verified by Kinnersley and Bonnor that the energy- 
momentum tensor (3.3) entails an outgoing flux of energy and momentum at 
infinity. They then argued heuristically that this must be compensated by 
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a loss of mass dm{s)/ds, and a recoil du^/ds of the "central particle". Let 
us verify this more rigourously by studying whether the total energy tensor 
T^^{x) can be everywhere conserved, including on C This will prove the con- 
sistency of the point-like limit. We shall consistently work in the framework 
of distribution theory, with T^"{x) considered as a linear functional acting on 
smooth, compact-support functions in M"^ : (T'^'^, cp) = J d^x T^'^{x) ^{x). For 
instance, {T^^^^ip) = J ds m{s) ip{z{s))^ while the distribution Tj^^^x) 

is uniquely defined from the function (3.3), in spite of its singular behaviour 
on the central worldline £, because Tj-^^ is locally integrable [0{p^^) singula- 
rity only as p ^ 0]. The distributional derivative of the matter part (3.2) is 
easily verified to give 

du TlZ^ix) =+Jds-^^ {m{s) u^) 6,{x - z{s)) . (3.4) 

The computation of T^^^^ is more tricky. It cannot be directly done on 
the function (3.3) because d^, T^p^(x) has a non locally integrable singularity, 
0(p~^), along C One way to do it is to apply the distributional definition, 
(^(p)' ~^>' Stokes formula, with equation (2.12), when integrating 

by parts on a domain Pr{x) > e before letting e — 0. A more convenient 
technique, which is quite standard in distribution theory for functions ha- 
ving power-law singularities (see e.g. @]), is to introduce a complex para- 
meter B and to work with the analytic continuation of integrals of ordinary 
functions. For instance, we can multiply equation (3.3) by p^ and compute 
dy{p^ T^J^{x)) which is a locally integrable function for Re{B) > 0. The 
looked for distribution is then obtained by analytically continuing B to zero. 
Using the derivatives recalled in section 2, this gives after an easy calculation 

d^ip^ Ti;^) = B p^-' ^^^1^ F . (3.5) 

The problem is then reduced to evaluating the analytic continuation down 
to i? = of the four- dimensional integral 

/(5) = J d^x d^{p^ Tjp ^{x) = ^Jdsdp dilr, p^~^ eis,n) F ifi{x) , 

(3.6) 

where we have used retarded polar coordinates and equation (2.11). The B 
factor in front shows that only an arbitrarily small neighbourhood of £, say 
p < pq, matters. Using the Taylor expansion of ip{x) gives ip[x{s, p,9,(p)] = 
ip[z{s)] + 0{p) when p — >• 0. As the analytic continuation of B /q" dp p^~^ p = 
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B Pq^^/{B + 1) is zero, we see that only the value of ip along C, ip[z{s)], 
contributes. The integration over the angles introduces the quantity 

/ dn^ e{s, n) P = eo{s) + eI{s) (3.7) 

with 

eo{s) = ^Jdn^e{s,n), (3.8) 

e^(s) = — fdQ^ e(s, n) , (3.9) 
An J 

while the integration over p gives, modulo terms that tend to zero when B 
is analytically continued to zero, B Jq° dp p^^^ = Po = e^''^'"', which tends 
to one when B is continued down to zero, independently of po- Finally, this 
gives 

/(O) = {d, T(;^, ip) = J ds{eo{s) + e^,{s)) ^[z{s)] , (3.10) 
which means that we have proven the distributional result 

du Tj^(x) = j ds{eo{s) + et{s)) 5,{x - z{s)) . (3.11) 

Combining (3.11) with (3.4), we conclude that the total energy-momentum 
tensor (3.1) will be everywhere conserved in the sense of distribution theory, 
{d^ T^^ (^) = 0, if and only if 

(m(s) u^{s)) + e^{s) + £^(s) = . (3.12) 

ds 

Noting that Si = from its definition (3.9), equation (3.12) is equivalent 
to 

m = —Eq , (3.13a) 

mil'' = -e^ . (3.13b) 

These are the results of obtained here without relying on a heuristic 
energy- momentum balance argument. 

4 Linearized gravitational waves emitted by 
photon rockets 

The conserved energy-momentum tensor (3.1) can be consistently taken as 
distributional source term for linearized gravity. The linearized Einstein 
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equations for g^y = rj^^ + h^,, (or the field equations of the massless spin-2 
field h^u) read 

□ h^, + d,, h - d,^ K - d,^ hi = -1671 G {t,, - ^ r^,, t) , (4.1) 

where □ = 77"^ S^/j, h = h^, T = , indices being raised and lowered by 
the Minkowski metric. To make contact with usual perturbation techniques 
in general relativity let us start by looking for solutions of equation (4.1) in 
harmonic gauge 

/i?."^ h^^""'^ = . (4.2) 

This yields as usual 

□ = -levr G (t,, - \ V,. • (4.3) 

To simplify the discussion we assume that the photon rocket started (smooth- 
ly) its activity only at some finite time in the past, i.e. e{s, n^) = for s < sq. 
This implies that m(s) was constant and C straight before sq (see equations 
(3.13)), and that the non compact-support part of T^^{x), namely T^^^{x), 
equation (3.3), has its support entirely contained in the forward light cone 
with vertex z^{sq). These properties ensure that the retarded solution of 
(4.3), i.e. the convolution of the retarded Green's function, 

Gnix'') ^ d{x') 5{^,, x^) ^ - . □ G'^(^) ^ 

2tt An |x| 

(4.4) 

with the right-hand side of (4.3), is mathematically well- defined]^, and can 
be characterized as being the only solution of (4.3) containing no incoming 
radiation. We write it as 

h%^^\x) = -16vr G n-^' (t,^ - \ t) , (4.5) 

where {O]^^ S){x) = J d'^y Gr{x — y) S{y) denotes the retarded potential 
operator. 

The linear split (3.1) of T^^, leads to a corresponding linear split of /il^Jf^^ 

say 

h^rHx) = h^;:\x) + h^^}{x) . (4.6) 

■^In particular, it is given by a compact-support three-dimensional integral; the "bad" 
~ 1 /r^ behaviour of the source term at future null infinity causing no convergence problem 
in the retarded potential integral. 
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The "matter" part /i^™-* is obtained in a standard way from its source (3.2) 
using the covariant form of Gr{x^ — z^{s)) oc 5[{x^ — z^{s)){x^ — z^{s))] 
and the formula 5[<y3(s)] = ^{s — Sn) /\^'{sn)\ where n labels the roots s„ 

n 

of (^(s) = 0: 

/.ff(x) = 2G ^mu,u. + U^ (4^7) 

Here, as in section 2, the suffix R indicates that all s-dependent quantities 
must be evaluated at the retarded point sr{x) (Lienard-Wiechert potential). 

The "photon" part /ijfj is given by the following retarded potential inte- 
gral 

h^^ix) = -AG U,- [^£(£i^^^^ (4 g) 

which is well-defined as the integral of a function (rather than a distribu- 
tion) because of the local integrability of The integral (4.8) can be 
analytically reduced to a one-dimensional integral over the proper time s (on 
the domain s < sr{x)) by using the formula (2.24) of Ref. []10|Q- However, 
a more convenient way of studying h'^^} is, following the study of the m"^ 
term in Ref. [|T^, which has the same structure as h'^^} (for a straight world- 
line), to introduce a simplifying gauge transformation. The role of this gauge 
transformation (or linearized coordinate transformation) is to transform the 
source oc k^j (? in equation (4.8) into a source proportional to which 
turns out to be much simpler to deal withQ. However, in so doing we have to 
tackle non locally integrable terms ~ . This situation is technically very 
similar to what happened in section 3 where we needed to evaluate the dis- 
tributional derivative of Tji^^ ~ p~^, i.e. to make sense of integrals involving 
d\ T^p^ ~ p~^. Like in section 3, a technically convenient, and mathema- 
tically rigourous, way of dealing with such terms is to introduce a complex 
parameter B and to work with the analytic continuation of certain integrals. 

Given a function (t(s, n'^) which reduces to an angle-independent constant 
(To for s < So (i.e. when £:(s,n'^) vanishes) and which will be related to m(s) 

'^To generalize the result of which assumes a straight worldline L to the present 
case of a curved worldline, one needs (beyond multipole-expanding e(rt)) to represent the 
source term in (4.8) as a one-dimensional integral of elementary contributions with support 
on light cones with vertex on L. 

^Another utility of that gauge transformation is to deal away with the terms of order 
(In p)l p in /i^'jf'^'' at future null infinity. As is well-known these terms, due to the slow 
fall off of T^i/, arise in harmonic coordinates but can be gauged away by introducing 
some suitable, '"radiative" coordinates. 
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and e{s,n'^) below, we define the gauge transformation ^^(x) as follows. We 
first consider the retarded potential integral 

(x) = a-' (^1 a{sn{x),n) k,^ , (4.9a) 

which depends on the complex parameter B, and is convergent when — 1 < 
Re{B) < 0. Then we analytically continue ^^(x) in 5. Generalizing standard 
arguments (see, e.g., Refs 0, ^) ^^(x) is easily seen to admit a continuation 
as a meromorphic function over the whole complex B plane. Denoting by 
FPb=q ("Finite Part at B = 0") the operation of taking the constant term 
(zeroth power of B) in the Laurent expansion of ^^(x) around = 0, we 
then define 

Ux) ^ FPb=o (ef (^)) • (4.9b) 

As is locally integrable the analytic continuation factor is not needed 
in (4.9a) to deal with the neighbourhood of C, but is useful to deal with 
the slow fall-off (~ p~^) at past null infinity due to the fact that a becomes 
a non-zero constant in the past. The retarded integral (4.9a) is convergent 
when —1 < Re{B) < 0, and is easily seen to have a simple pole at i? = 
when analytically continued upwards in Re{B). The effect of this pole is to 
introduce a term proportional to do mJ^ Pq / B in and therefore, using the 
Laurent expansion, B^^ p^ = B^^ exp{B In po) = + In po + 0{B), a 
term proportional to ctq m° In po in [Here, po denotes the spatial distance 
between the field point x'^ and the incoming straight worldline with four- 
velocity This logarithmic term disappears when taking x-derivatives, 
which corresponds to the fact that the 5-dependent integrals considered 
below will have a better convergence at past null infinity, and generate no 
poles at i? = 0. 

Thanks to the properties of analytically continued integrals we can com- 
mute x-derivatives with the integral operator O^^. [When dealing with an 
integral such as (4.9a), the integrand of which has power-law singularities 
both on C and at past null infinity, it is convenient to split the spatial do- 
main of integration in two pieces, one of them being compact and centered 
on C] see, e.g., section 3 of 0| for a similar situation.] This yields 

e = I [Mp""-' ^ k^)] , (4-10) 
where, using the results of section 2, 

^/^(pI^^ 0^(5/2, n") K) = -p|"^ Da k^ 
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+PR~^ Y^i^nu + {B-2) u^ + {B- 3) n^, n^} 

) K\ . (4.11) 

The quantity Da in (4.11) denotes the action of a certain first-order diffe- 
rential operator on a{s,n"): 

Da{s, n") = + [{nil) K - ^ + (3 - B){nu)a . (4.12) 

Here we consider (t{s, rf) as some exphcit function of s and the four compo- 
nents of n", e.g. the one obtained by expanding the angular dependence of 
a in cartesian symmetric-trace-free (STF) polynomials in n° [equivalent to 
a spherical harmonic expansion in Yf^{9, (p)] 

oo 

a{s, n") = J2 ^a,a,...ae {s) n^i----^ , (4.13) 

l=Q 

where n"^ " "^ denotes the STF projection (within the spatial hyperplane or- 
thogonal to w") of . . .n"^ For instance, n^"' = n^'n''-l A^"". Without 
loss of generality, we can also consider that the s-dependent coefficients in 
(4.13) are purely spatial (-u^^ o"q,^...q,^ = . . . = u"^ (Jai...ai = 0) and STF, 
in which case we can replace n"^ ' "* by n"^ . . .n°'^ in (4.13) so that the s- 
derivative in (4.12) acts only on Ua^ . Q^ (s). Evidently, the quantities Da 
and da / dna{A.a^ — Ua entering (4.11) are defined independently of the 
multipole expansion (4.13), but an alternative definition would imply choos- 
ing some explicit parametrization of n" in terms of two polar angles 9^ ip (i.e. 
choosing a triad of vectors orthogonal to -u") which we prefer to shun doing. 

We see from equation (4.11) that the gauge transformation h'^^ = h^^i, + 
dfj, + dy generates terms proportional to k^,/ which can cancel the 
original "photon" piece (4.8) and replace it by the retarded potential of source 
terms oc p~^. To perform this transformation explicitly, it is convenient to 
consider that the photon flux function e[s, n") is, similarly to equation (4.13), 
decomposed in angular multipoles 

oo 

£(5,n")= 1:^(5, n") , (4.14) 

where ee{n") is a linear combination of f^;"!---"*. The first two terms in the 
multipole expansion (4.14), namely the monopole Eq and dipole parts £i, are 
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fully determined by the integrals (3.8) and (3.9). More precisely eo (defined 
by equation (4.14) as the monopole piece of e{s,n)) is identical with eo{s) 
defined by equation (3.8), while the dipole piece of e{s,n) reads 

£i(s,n") =3£^(s) (4.15) 

in terms of the definition (3.9). 

Let us first consider the special case (that we shall the "Kinnersley case" ) 
where the photon flux e{s, n) contains only a monopole and a dipole piece: 

£^(s,n") =£o(s) + 3£^(s) n„ . (4.16) 

In view of the balance equations (3.13), this can be rewritten in terms of 
"mechanical" quantities 

£^(s, n") = -[m(s) + 3 m{s) {ii n)] . (4.17) 

Let us correspondingly consider the special case where a is purely monopo- 
lar: a{s,n) = a^{s). From equations (4.10)-(4.12) (which simplify when a 
depends only on s) we find for the gauge transformation 

d, + 5. = n^^j-p^-^ [a^ + (3 - S) (nu)] k, K 



]^ 2 

— H — (n^ Uu + rii, U/j) + {B — 3) 



(4.18) 



where we have decomposed the functions of in irreducible multipoles of 
order 0, 1 and 2 {n^''' = 7%^" rf -\ A^""). 

Comparing (4.18) with the result of inserting equation (4.17) into equa- 
tion (4.8) we are led to defining 



a"^ (s) = +A G m{s) . (4.19) 

For this choice we see easily that the first term on the right-hand side of 
(4.18) is analytic at i? = (no pole; thanks to the vanishing of and 
a^{nu) in the past) and that its value at = is 

which is precisely the opposite of h'^^{x)^ defined in equation (4.8), in the 
Kinnersley case. 
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We are then left with evaluating three types of terms: "monopolar" terms 
of the form p^~^ and "dipolar" and "quadrupolar" terms of 

the form □]j^(p^~^ (^(s) Ua^^^ at) with £ = 1 or 2. The monopolar term 

-1 (B R-3 



Pr'^'M AM.(^K)j (4.20) 

appearing in the second square bracket on the right-hand side of (4.18) is 
easily evaluated, when B is analytically continued to zero, by the same rea- 
soning used in section 3 above to deal with the integral of the product of the 
right-hand side of equation (3.5) with a smooth function [here that smooth 
function is G r{x ~ x')]. Denoting ACb=o the operation of taking the analytic 
continuation at i? = 0, we deduce from the results of section 3 

ACb=o {B p^-' ^{s)} = A'k j ds ^{s) 5^{x - z{s)) , (4.21) 

in the sense of distribution theory [where the B factor in front ensures the 
absence of pole at B = 0]. Convoluting the Green's function (4.4) with 
(4.21), we find that the analytic continuation of the monopolar term (4.20) 
yields 



1 a^jsR) A^,{sr) ^ AG ( m{s) A,,/ 
3 PR 3 I p ^ 



(4.22) 

R 



The dipolar and quadrupolar terms (last two terms on the right-hand 
side of equation (4.18)) can be evaluated by using explicit formulas derived 
in Ref. |^. There, indeed, it was shown that, in the case of a straight central 
worldline Cq, the analytic continuation at S = of the retarded potential 
generated by r^~^ H(t — r) ni^,„it was expressible, when 3 < < £ + 2, as 
a finite sum of derivatives of H{t — r) divided by powers of r. In particular, 
when k = 3, £ > 1 and H(t) is the one- dimensional delta function 6(t — a), 
equation (4.24) of Ref. reads 

ACb=o {On\r''-' 5{t ~ r - a) = 

Here, the central straight worldline £o is taken as time axis, r denotes the 
spatial distance between the field point {x^) = (t,x*) and Cq, ra* = x'^/r, and 
0" is a constant, parametrizing a specific point {a, 0) on Cq. 
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Equation (4.23) can be generalized to the case, of interest for the present 
paper, of a generic curved worldhne C by noticing that if Co is taken to be the 
tangent to C at z^{a) then [S{t — r — o-)]co — [S{sr{x) — a)]c as distributions 
in R^. [Indeed, both distributions have the future hght cone of z^{a) as 
support, and using the fact that the volume element in retarded coordinates 
(2.11) does not depend on the acceleration -u^ one sees immediately that 
/ d^x 6{t — r — a) (p{x) = J d*x 5{sr{x) —a) ^{x)]. Having generalized (4.23) 
to an elementary source with support on a light cone centered on £, we can 
pass to a general retarded-time-dependence of the source by integrating over 
the parameter cr using the decomposition H{sr{x)) — J da H{a) S{sR{x)—a). 
This yields 

ACb=o {a-n' (pr' H{sn{x)) n,,...,,)} = 

L_ ^('^^(^)) fi ■(£>!) (4 24) 

It is remarkable that, though the left-hand side, evahiatcd at the field point 
xf^, is an integral over the past light cone of whose integrand depends on 
the full past history of the source, i.e. on the value of the function H(s) over 
the interval —oo < s < sr{x), the final result on the right-hand side depends 
only on the value of H{s) at the retarded point sr{x). This simplification 
does not hold for integrands proportional to p^^, as they appeared in the 
original expression (4.8) before introducing the gauge transformation (4.9). 

Using the result (4.24) we can read off equation (4.18) the contributions 
of the dipolar and quadrupolar terms in the last bracket. Namely: 



PR 



+ 2 ("-M ^'^ + "''^ + 2 ^^'^ 



2G ^^^^^ [n^ + n^u^ + V] ■ (4-25) 

PR 

Gathering our results, we conclude that the gauge-transform of the original 
harmonic linearized field (4.6), 

Kr^i^) = + 9m + = h^;:\x) + /^^V) , (4.26) 

with 

htl (^) = ^ (^) + d,i. + , (4.27) 

is given, in the Kinnersley case (4.16), (4.17), by adding the matter contribu- 
tion (4.7) and the gauge-transformed photon one obtained by adding (4.22) 
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and (4.25): 



K 



2G 



PR 



R 



(4.2^ 



Replacing the definitions n^^, = 
eliminating in favour of A;^ 

as 



n,. 



A 



Uf, + Hi,, 



= Vf^u + and 
the result (4.28) can be rewritten 



- i A 

3 ^Mi'' 



K 



2G 



m[SR) 
PR 



2 Ui 



(4.29) 



We recognize in the first two terms of the (gauge-transformed) photon contri- 
bution (4.29) the opposite of the matter contribution (4.7). We have therefore 
proven by an explicit computation that, after a suitable gauge transforma- 
tion (or linearized coordinate transformation) (4.9), there were terms in the 
gravitational field generated by the photon energy-momentum tensor which 
cancelled the usual, direct matter terms (4.7) to leave as final combined 
gravitational field in the special Kinnersley case (4.16), (4.17) the net result 



K 



L(har) 
"'fii/ 



K 



K 



2G 



m[SR) 
PR 



ky 



(4.30) 



Because of the specific algebraic structure oc k^ k^ of the last right-hand side 
of equation (4.30), the gravitational field is easily verified not to contain 
any physical 1/r gravitational waves at (future null) infinity. This can be 
seen either by taking the usual transverse-traceless (TT) projection of the 
spatial components in some fixed reference frame [using fcj = nj + 0(l/r)], 
or by checking, from the derivatives given in section 2, that, at the 1/p level, 
dai3 OC ka kp k^ ky which implies the vanishing of the linearized Riemann 
tensor at the 1/p level. [A more complete calculation, first performed by 
Kinnersley |jl|, shows that the linearized Weyl tensor of (4.30) is of order 
at infinity, while its linearized Ricci tensor, proportional to T^y, is of order 

p-\ 

The non radiative result (4.30) has been derived in the special case where 
the anisotropy of the photon fiux e(s,n) is purely dipolar, equations (4.16), 
(4.17). Let us briefiy discuss the generic case where e{s,n) contains more 
general anisotropies. Using the notation (4.14), we can always write 



n) 



n) + 



e=2 



,n) 



(4.31) 



16 



keeping unchanged the definition (4.16), and the result (4.17) which followed 
from the general balance equations (3.13). Let us then define a{s,n) as the 
unique solution of the differential equation 

Do a{s, n) = -4 G e{s, n) (4.32) 

which reduces to the constant 4 G m(so) when s < Sq (i.e. before £{s, n) turns 
on). Here Dq denotes the value at S = of the differential operator (4.12). 
The (local) existence and uniqueness of a follow from the fact that, when 
written explicitely in terms of s and two angular coordinates {6 a) = {0, <^) 
on the sphere, Dq = d/ds + v^{s, 6b) d/dOA is a first-order linear operator. 
Given initial data for C in the past, and the monopole and dipole parts, Sq 
and £i, for all values of s one can solve uniquely for the time evohition of C 
(via equations (3.13)) and then for a{s^n). The linearity of equation (4.32) 
means that its solution reads a = + a'^^'^, where is the "Kinnersley" 
piece (4.19) and where a'^^ denotes the solution of Dq cr""^ = — 4 G J2 ^i- 

From the explicit expression (4.12), one checks that, when expanding a'^^ in 
orbital multipoles oc nai. the £ = and i = \ pieces of Dq g^^^ are 
entirely expressible in terms of crj^ and (without the presence of a term 
oc a^^ n^,). As the source term for Dq a^^'^ has only multipoles of order 
> 2, we conclude that Gq"^ = af^ = 0, i.e. 

oo 

a{s,n) = AG m{s) + ^ai{s,n) , (4.33) 

where (72 + cxs + ■ ■ ■ is driven by £2 + ^3 + ■ • ■• 

Using as above a to define the gauge transformation (4.9), we find from 
equation (4.11) that the total gauge-transformed gravitational field reads 

where [h' is the non-radiative field (4.30) and 



h' 



(4.34) 
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^(p)rad 



1 da'^ 



B 



A 



B-2 



{rifj, Uu + Uy Uf,) + {B -3) n 



2 dn^ 

From formula (4.24), which can be compactly stated as 



AC 



B=0 



□ 



-1 
R 



pI-' H{sr{x) 



(4.35) 



(4.36) 



where A„ is the Laplacian on the unit sphere, and A„ its inverse acting 



pR U^y{sR,n) 



where 



on zero-mean functions on the sphere, we get h^j'^}^'^'^ 
values of for the same value of the proper time s. This computation has 



n) is some explicit functional of (j'''^*^(s, n), which depends only on the 



see 



been done explicitly in the case of a straight worldline C in Ref. [|10[ 
equations (2.20) there. In particular, the TT projection of equation (2.20c) 
there [in which only the last term oc survives] show that the 

1/r gravitational wave amplitude radiated at infinity is non zero as soon as 



cr 



rad / 



n 



contains a non-zero multipole of order l>2. 



The conclusion is that, as soon as the photon flux function e(s, n) contains 
an anisotropy of multipole order £ > 2, it will drive such terms in a'^^{s,n) 
which will, through equations (4.35), (4.36), radiate genuine 1/r gravitational 
waves at future null infinity^. 



5 Conclusion 

After having shown the mathematical and physical (Vx G M'*, dp T^'^{x) = 0) 
consistency of the definition of massive point-like photon rockets in special 
relativity, we have studied the retarded linearized gravitational field that 
they generate. The crucial point is that the source of the gravitational ra- 
diation is the sum of the stress-energy tensor of the matter and of that of the 
outgoing flux of photons. In fact, the basic source of time-dependence in the 
system is the photon flux function e{s,n). The acceleration of the matter 

^Note that in the Umiting case of an infinitely massive photon rocket, £ becomes 
straight and a^'^^{s, n) ~ —4 G ds' ^ ££(s', n), so that we can directly use the results 

°° t>2 

of to relate the multipoles of the gravitational wave at infinity to the multipoles of 
order ^ > 2 in e(s, n). 



18 



is a derived quantity, computable from the monopolar and dipolar pieces 
in the spherical harmonics expansion of e{s,n): equations (3.13). We have 
proven by an explicit calculation of the retarded potential generated by the 
total stress-energy tensor T'^" = T^^^ + Tj^^ that the amount of gravitational 
radiation at infinity is driven by the anisotropics of order £ > 2 in the photon 
flux e{s,n). On the other hand, the coherent addition of the gravitational 
waves emitted by the monopolar and dipolar pieces of the photon energy- 
tensor, and of the gravitational waves emitted by the accelerated (because of 
the recoil) matter, exhibits a cancellation leaving a time-dependent, but non- 
radiative gravitational field (4.30). This cancellation is not too surprising if 
wc consider, as recalled above, that the photon fiux function e{s,n) is the 
basic source of time-dependence, and therefore gravitational radiation, in the 
system. In view of the spin 2 of the graviton, it makes sense that only the 
anisotropics of order ^ > 2 in e{s, n) can radiate gravitational waves. 

We have presented the calculation of the emission of gravitational radia- 
tion in gory details to make it clear that there is absolutely no incompatibility 
between the non-radiative character of the gravitational field generated by a 
purely dipolar photon flux anisotropy and the perturbation techniques cur- 
rently used to compute the generation of gravitational radiation in general 
relativity. Some of the technicalities of our treatment (notably the use of 
analytic continuation techniques) have been forced upon us by the necessity 
to deal in a mathematically rigourous way with the retarded potentials gene- 
rated by peculiar singular sources oc p]^ (which become inflnite on a world- 
hne and fall off rather slowly at infinity). A careless treatment could miss 
the delicate compensation between the gravitational radiation emitted by 
the Dirac-distributed matter energy tensor, and that emitted by some of the 
singular terms oc appearing as effective sources in the gauge-transformed 
photon-generated metric /i^^^ + d^j, + di, 

Appendix 

To make contact with the literature let us remark that, as a by product of 
the calculations above, we have shown that the linearized metric 

= V,. + C = v,. + '2G K , (A.i) 

Pr 

satisfies everywhere in ^ (in the sense of distribution theory) the linearized 
Einstein equations 

R'il - \ V^.. = SttG 77.A Ti , (A.2) 
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with given by equations (3.1)-(3.3). 



Remarking now that (A.l) is of the Kerr-Schild form [12|, and remem- 



bering that Giirses and Giirsey |TB[ have shown that all the nonlinear terms 
cancel when Kerr-Schild metrics are inserted into the Einstein equations writ- 
ten in terms of a mixed energy- momentum tensor T^, we see that we have 
rederived (in a roundabout way) the result of Kinnersley namely that 
the metric (A.l) is, outside the worldline £, an exact solution of Einstein 
equations with source Ty\n = (4vr)~^ p^^ e^{sji,n) . Note that 



m 



the exact, nonlinear framework of Einstein's equations, it makes no sense to 
write down 5-like source terms, and therefore to say that the exact Kinners- 
ley metric describes a point-like photon rocket. This is why, in the text, we 
consistently considered only linearized gravity. 



20 



References 

[1] Kinnersley W 1969 Phys. Rev. 186 1335. 

[2] Kramer D, Stephani H, MacCallum M and Hcrlt E 1980 Exact solutions 
of Einstein's field equations (Cambridge: Cambridge University Press). 

[3] Bonnor W B 1994 Class. Quantum Grav. 11 2007. 

[4] Schwartz L 1966 Theorie des distributions (Paris: Hermann). 

[5] Hornby A S 1974 Oxford's advanced learner's dictionary of current En- 
glish (Oxford: Oxford University Press). 

[6] Blanchet L and Damour T 1986 Phil. Trans. R. Soc. Lond. A 320 379. 

[7] Blanchet L and Damour T 1988 Phys. Rev. D 37 1410. 

[8] Blanchet L and Damour T 1989 Ann. Inst. H. Poincare 50 377. 

[9] Damour T and Iyer B R 1991 Ann. Inst. H. Poincare 54 115. 
[10] Blanchet L and Damour T 1992 Phys. Rev. D 46 4304. 
[11] Chown M 1994 New Scientist 143 No 1943 ( 17 September 1994 ) p.l5. 
[12] Kerr R P and Schild A 1965 Proc. Symp. Appl. Math. 17 199. 
[13] Giirses M and Giirsey F 1975 J. Math. Phys. 16 2385. 



21 



